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Abstract. In interacting multi-component dark matter (DM) models, if the DM components 
are nearly degenerate in mass and the interactions between them are strong enough, the 
relatively heavy DM components can be converted into lighter ones at late time after the thermal 
decoupling. Consequently, the relic density of the lightest DM component can be considerably 
enhanced at late time. This may contribute to an alternative source of boost factor required to 
explain the positron and electron excesses reported by the recent DM indirect search experiments 
such as PAMELA, Fermi-LAT and HESS etc.. 



In the recent years, a number of experiments such as PAMELA [T], ATIC [5], Fermi-LAT [3] 
and HESS [1] etc. have reported excesses in the high energy spectrum of cosmic-ray positrons 
and electrons over the backgrounds estimated from the standard astrophysics, which may be 
interpreted as indirect signals of the annihilation or decay of dark matter (DM) in the Galactic 
halo. If the DM particles are thermal relics such as the weakly interacting massive particles 
(WIMPs), the thermally averaged product of their annihilation cross section with the relative 
velocity at the time of thermal freeze out is typically {(Tv)p ~ 3 x 10"^^cm^s~^. The positron 
or electron flux produced by the DM annihilation can be parametrized by 



where po is the smooth local halo DM energy density estimated from astrophysics, A^e is the 
averaged electron number produced per DM annihilation which depends on DM models and 
parameters in the models for the propagation of cosmic ray particles, and is the mass of 
the DM particle. The boost factor B is defined as B = {p / po)^ {av) / {av) p with p the true local 
DM density and {av) the DM annihilation cross section multiplied by the relative velocity and 
averaged over the DM velocity distribution today. Both the PAMELA and Fermi-LAT results 
indicate that a large boost factor is needed OE]. For a typical DM mass of ~1(1.6) TeV the 
required boost factor is i? ~ 500(1000) for DM annihilating directly into p^ and p fixed at 
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Po = 0.3 GeV • cm 
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A large boost factor may arise from the non-uniformity of the DM distribution in the halo. 
The N-body simulations show, however, that the local clumps of dark matter density are unlikely 
to contribute to a large enough p/po [ZllE]- An other possibility of enhancing the boost factor 
is that the DM annihilation cross section may be velocity-dependent which grows at lower 
velocities. The DM annihilation cross section today may be much larger than that at the time 
of thermal freeze out, and thus is not constrained by the DM relic density. Some enhancement 
mechanisms have been proposed along this line, such as the Sommerfeld enhancement [5Vll7|. 
In some non-thermal DM scenarios, the number density of the DM particle can be enhanced by 
the out of equilibrium decay of some heavier unstable particles if the DM particle is among the 
decay products of the decaying particle [TS1[T2]. The decay of the unstable particle must take 
place at very late time. Otherwise the DM particles with the enhanced number density will 
annihilate into the Standard Model (SM) particles again, which washes out the enhancement 
effect. 

In this talk, we discuss an alternative origin of the boost factor arising from the late time dark 
matter conversion processes, which requires neither the velocity-dependent annihilation cross 
section nor the decay of unstable particles [20]. We show that in the scenarios of interacting 
multi-component DM, the interactions among the DM components may convert the heavier 
DM components into the lighter ones, which is not sensitive to the details of the conversion 
interactions. If the interactions are strong enough and the DM components are nearly degenerate 
in mass, the conversion can enhance the number density of the lighter DM components at late 
time after the thermal decoupling. Eventually, the whole DM today in the Universe may consist 
of only the lightest DM component with enhanced number density, which leads to a large boost 
factor. The scenarios of multi-component DM have been discuss previously in Refs. |21H30j . 
Note however that the models with simply mixed non-interacting multi-component DM cannot 
generate large boost factors. 

Let us consider a generic model in which the whole cold DM contains components 
Xi {i = 1,...,N), with masses rrii and internal degrees of freedom gi respectively. The DM 
components are labeled such that rrii < irij for i < j, thus xi is the lightest DM particle. 
We are interested in the case that Xi nearly degenerate in mass, namely the relative mass 
differences between Xi &iid xi satisfy £i = {mi — mi) /mi <^ 1. In this case, we shall show 
that the interactions between the DM components lead to the DM conversion. The thermal 
evolution of the DM number density normalized to the entropy density Yi = rii/s with respect 
to the rescaled temperature x = mi/T is govern by the following Boltzmann equation 



dYi{x) _ A 
dx x^ 



{aiv){Y^ - y^) - Y.{a,,v){Y^ - r^Y^ 



(2) 



where A = xs/H[T) is a combination of x, the entropy density s and the Hubble parameter H(T) 
as a function of temperature T. Yi^q ~ (^j/s) [mjT/ (27r)]^/^ exp(— e^x) is the equilibrium number 
density normalized to entropy density for non-relativistic particles, {uiv) are the thermally 
averaged cross sections multiplied by the DM relative velocity for the process XiXi ~^ XX' with 
XX' standing for the light SM particles which are in thermal equilibrium, and {(Jijv) are the 
ones for the DM conversion process XiXi ~^ XjXj- The quantity 



Yjeq(x) 



= J^Mn = ( 7^ ) ( ) exp[-(e, - ej)x] (3) 



is the ratio between the two equilibrium number density functions for components i and j. In 
Eq. ^ we have assumed kinetic equilibrium. The first term in the r.h.s. of Eq.(l2|) describes the 



change of number density of Xi due to the annihilation into the SM particles, and the second 
term describes the change due to the conversion to other DM components. 

In the case that the cross section of the conversion process {dijv) is large enough, the DM 
particle Xi can be kept in thermal equilibrium with Xj for ^ long time after both Xi ^^'i Xj 
have decoupled from the thermal equilibrium with the SM particles. In this case, the number 
densities of Xi,j satisfy a simple relation 



Uj{x). (4) 



Even when Xi is in equilibrium with Xj the ratio of the number density Yi{x) /Yj (x) can be quite 
different from unity and can vary with temperature. For instance, if gi ^ gj and < (ei—ej) <^ 1, 
from Eq. ([3D and (jH) one obtains Yi{x) S> Yj{x) at the early time when (e, — ej)x <^ 1. However, 
at the late time when (e^ — ej)x ^ 1, one gets Yi{x) <C Yj{x), which is simply due to the 
Boltzmann suppression factor exp[— (e^ — ej)x] in the expression of rij. Thus the heavier particles 
can be gradually converted into lighter ones through this temperature-dependent equilibrium 
between Xi and Xj- 

An interesting limit to consider is that the rates of DM conversion are large compared with 
that of the individual DM annihilation into the SM particles, i.e. {(Jijv) ^ {civ). In this 
limit, after both the DM components have decoupled from the thermal equilibrium with the SM 
particles, which take place at a typical temperature x = x^ec ~ 25, the strong interactions of 
conversion will maintain an equilibrium between Xi and Xj for a long time until the rate of the 
conversion cannot compete with the expansion rate of the Universe. Making use of Eq. ([3]), the 
evolution of the total density Y{x) = Yl,i=i can be written as 

f = -^(-e//.)r-y4), (5) 

where {aeffv) is the effective thermally averaged product of DM annihilation cross section and 
the relative velocity which can be written as 

/ \ E^i Widfjl + Sj)^ exp(-2£^x) 

^eff 

where Wi = {aiv) / {(Jiv) is the annihilation cross section relative to that of the lightest one. The 
total equilibrium number density can be written as 

Yeq = 5^ Yieq{x) ^ geff (j^J exp(-x), (7) 

with an effective degrees of freedom g^^ff = ^ • gi{l + Cj)^/^ exp(— ejx). Note that the conversion 
terms do not show up explicitly in Eq. ([5]). Through the conversion processes XiXi ~^ XjXj the 
slightly heavier components will be converted into the lighter ones, because the factor rij[x) is 
proportional to exp[— (mi — mj)/T] which suppresses the density of the heavier components at 
lower temperature. If the conversion cross section is large enough, most of the DM components 
will be converted into the lightest xi before the interaction of conversion decouples, which may 
result in a large enhancement of the relic density of xi and leads to a large boost factor. 

As an example, let us consider a generic DM model with only two components. For relatively 
large conversion cross section u = {(^21^) / {ctiv) > 1, The effective total cross section is given by 
{cjeffv) = {aiv)[l + w^f^ exp(— 2ex)]/[l + 51 exp(— ex)]^], where w = W2, g = 92/91 and e = £2- 



Because of the x-dependence in {a^ffv), the thermal evolution of Y{x) differs significantly from 
that of the standard WIMP. In the case that X2 has large degrees of freedom but a small 
annihilation cross section, namely 5 ^ 1, w <C 1 and wg"^ <C 1, the thermal evolution of the 
total density Y can be simplified. The thermal evolution of the total number density can be 
roughly divided into four stages: i) At high temperature region where 3 < x <C Xdec, both 
the DM components are in thermal equilibrium with the SM particles. Yi{x) must closely track 
Yieq{x) which decrease exponentially as x increases. However, since g ^ 1 and e <C 1, the number 
density of X2 is much higher than that of xi, i-e. Y2{x) » Yi{x). ii) When the temperature goes 
down and x is close to the decoupling point x^ec^ both the DM components start to decouple 
from the thermal equilibrium. In the region x^^ec ^ 2; ^ 1/e, {a^ffv) is nearly a constant and 
(aeffv) ~ {aiv)/{l + g)'^ <^ (fif), the total density Y{x) behaves just like that of an ordinary 
WIMP which converges quickly to Y{x) ~ Xdec/i^{(^iv))- hi) As x continues growing, the 
suppression factor exp(— ex) in {a^ffv) becomes relevant. The value of {a(,ffv) grows rapidly 
especially after x reaches the point ex ~ which leads to the further reduction of Y(x). 

In this stage, although both xi,2 have decoupled from the thermal equilibrium with the SM 
particles. The strong conversion interaction X2X2 ^ XiXi maintains an equilibrium between the 
two DM components. According to Eq. the relative number density Y2{x)/Yi(x) decreases 
with X increasing, which corresponds to the conversion from the heavier DM component into 
the lighter one. At the point Xc = (1/e) Ing one has Y2{x) Yi{x). For the region x > x^ec and 
X is not close to Xc, because of Yeq{x) <^ Y{x) and g exp{—ex) » 1, the total number density 
can be analytically integrated out, and Y{x) in this region can be approximated by 



iv) When x becomes very large ex » C(l) , Weffv) quickly approaches (o"if), and becomes 
independent of x again. The evolution of Y{x) in this region can be obtained by a simple 
integration as it was done in the stage ii). The solution of Y{x) shows a second decoupling. 
Finally when the conversion rate cannot compete with the expansion rate of the Universe at 
some point xp corresponding to sY2{cr2iv)/H ~ 1, both Yi{x) and Y2{x) remain unchanged as 
relics. The whole DM can be dominated by xi if the conversion is efficient enough. 

By matching the analytic solutions of Y{x) in different regions near the points Xdec and Xc, 
and requiring that the final total relic density is equivalent to the observed ^cDMh'^ ~ 0.11, we 
obtain the following approximate expression of the boost factor 



As expected, the enhancement essentially comes from the conversion of the degrees of freedom. 
Thus the maximum enhancement is g^. The two terms in the r.h.s of the above equation 
correspond to the reduction of Y[x) during the late time conversion stages. For large enough g, 
the boost factor can be approximated by -B ~ 5^/(1 + ^Q^Xdec/^'^Q)- order to have a large 
boost factor, a small e ^ In g/{g'^x dec) is also required. As shown in Eq. ^ the boost factor is 
not sensitive to the exact values of the cross sections as long as the conditions w <^ 1 and u ^ 1 
are satisfied. 

We numerically calculate the thermal evolution of Yi(x) and the boost factor without using 
approximations for a generic two-component DM model. The results for w = 10~^, u = 10 and 
e = 2 X 10~^ is shown in Fig.[TJ The value of {a2v) is adjusted such that the final total DM relic 
abundance is always equal to the observed value ^cDAih'^- The mass of the light DM particle is 
set to mi = 1 TeV. For an illustration the ratio between the internal degrees of freedom is set to 
be large g = 60. From the figure, the four stages of the thermal evolution of Y{x) as well as the 
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crossing point can be clearly seen. The crossing point at x = Xc ~ 2 x 10"^ indicates the time 
when the number density of xi start to surpass that of X2 and eventually dominant the whole 
DM relic density. In this parameter set a large boost factor B ~ (aiv) / {av) p ~ 585 is obtained 
which is in a remarkable agreement with Eq. ([9]) with error less than ~ 5%. For a comparison, 
in Fig. [T]we also show the cases without conversions. 




Figure 1. Thermal evolution of the number densities Yi(x) (red solid) and Y2{x) (blue 
solid) with respect to x. The solid (dashed) curves correspond to the case with (without) DM 
conversions. The green dotted curve corresponds to the sum of Yi and Y2, for parameters g = 60, 
mi = ITeV, e = 2 X 10~^, w = 10~^ and n = 10 respectively. 

The whole DM in the universe necessarily contains multiple components, as the lightest active 
neutrino already contributes to a small fraction of the DM relic density. It is easy to construct 
models with more stable neutrinos or neutrinos with lifetime longer than that of the universe. 
For instance, in fourth generation models with right-handed neutrinos, extra stable neutrinos 
may be the keV scale sterlile neutrinos and the heavy Majorana neutrinos which are stable due 
to additional symmetries [31] . For models with multiple DM components, it is possible that 
there exists interactions among the DM components which may lead to the conversions among 
them. In this talk we consider a simple interacting two-component DM model by adding to the 
standard model (SM) with two SM gauge singlet fermionic DM particles Xi,2- The particles 
Xi,2 are charged under a local U{1) symmetry which is broken spontaneously by the vacuum 
expectation value (VEV) of a scalar field (p. The corresponding massive gauge boson is denoted 
by A which may cause the reaction X2X2 ^ XiXi- The stability of xi,2 is protected by two 
different global U{1) number symmetries. An SM gauge singlet pseudo-scalar rj is introduced 
as a messenger field which couples to both the dark sector and the SM sector. In order to have 
the leptophilic nature of DM annihilation, we also introduce an SM SU{2)l triplet field A with 
the SM quantum number (1,3, 1) and flavor contents A = (5"'""'", 5"'', (5"). The triplet carries the 
quantum number B — L=2 such that it can couple to the SM left-handed leptons II through 
Yukawa interactions i^AiL, but cannot couple to quarks directly. The VEV of the triplet has to 
be very small around eV scale, which is required by the smallness of the neutrino masses. As a 
consequence, the couplings between one triplet and two SM gauge bosons such as S'^^W^W^, 
S^W^Z^ and S^Z^Z^ are strongly suppressed as they are all proportional to the VEV of the 



triplet, which makes it difficult for the triplet to decay even indirectly into quarks through SM 
gauge bosons |32ft37j . If r] has a stronger coupling to A than that to the SM Higgs boson H 
and 4> then the annihilation products of the dark matter particles xi,2 will be mostly leptons. 

The Lagrangian of the model can be written as £ = Csm + The new interactions in Ci 
which are relevant to the DM annihilation and conversion are given by 



A 3 Xi{i^-m)x^ + {D^,<^)HD^'(|))-m^,^^' 

-(W + CV^) [Tr(AtA) + k{H^H) + C{<PU)\ , = 1, 2) (10) 



Note that (p and rj do not directly couple to the SM fermions. After the spontaneous symmetry 
breaking in V{cl)), the scalar (f> obtains a nonzero VEV {(f>) = v^l\f2 which generates the mass 
of the gauge boson niA = gAV4,- At the tree level, the three components of the triplet 5^"'", 5+ 
and are degenerate in mass, i.e. m^++ = = = uia. 

We assume that X2 has large internal degrees of freedom relative to that of xi, i-e., 92 ^ 91, 
which can be realized if X2 belongs to a multiplet of the product of some global nonabelian 
groups. For instance 92 = ^92 with §2 =16, 8, and 4 if it belongs to the spinor representation of 
a single group of 5*0(8), 50(6) and 50(4) respectively. When X2 belongs to a representation 
of the product of these groups, its internal degrees of freedom can be very large. 

At the early time when the temperature of the Universe is high enough, the triplet A can 
be kept in thermal equilibrium with SM particles through the SM gauge interactions. The DM 
particles Xi s-i'^ in thermal equilibrium by annihilating into the triplet through the intermediate 
particle r]. The annihilation X2X2 t]* ^ 6^^6^^,6^5'^, 6^6^* is an s-wave process which is the 
dominant contribution . The ratio of the two annihilation cross sections is w = {y2 / Vif' {91 / 92) ■ 
It is easy to get a very small w provided that y2 <^ yi and gi ^ 92- In order to have a large 
enough {aiv) S> {(Jv)f the product of the coupling constants yifi must be large enough, or the 
squared mass of r] is close to s. The cross section of the conversion process X2X2 A* ^ XiXi 
is suppress by gi / 92 and also the phase space factor \/l — 4m^ / s when s is close to 4m2 at the 
vary late time of the thermal evolution. However, the cross section can be greatly enhanced if 
ruA is close to a resonance when the relation s ~ is satisfied. In the numerical calculations, 
we find that for the following selected parameters: mi = ITeV, e = 1 x 10""^, gi = 1, 92 = 60, 
ruA = 500 GeV, = 1.5 TeV, tua = 2.02 TeV, yi = 3, ?/2 = 0.07, /x/mi = 3, and gA = 2.5, 
the following ratio of the cross section can be obtained 

i(; ~ 1 X 10~^, u ~ 0.5, and {aiv) / {av) f — 500. 

In this parameter set the relative mass difference between niA and 2m2 is around 1%. The 
corresponding boost factor is i? ~ 500, which is large enough to account for the PAMELA data 
for the dark matter mass around TeV. 

In summary. We have considered an alternative mechanism for obtaining boost factors 
from DM conversions which does not require the velocity-dependent annihilation cross section 
or the decay of unstable particles. We have shown that if the whole DM is composed of 
multiple components, the relic density of each DM component may not necessarily be inversely 
proportional to its own annihilation cross section. We demonstrate the possibility that the 
number density of the lightest DM component can get enhanced in late time through DM 
conversation processes, and finally dominates the whole relic abundance, which corresponds to 
a boost factor needed to explain the excesses in cosmic-ray positron and electrons reported by 
the recent experiments. 
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